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We consider the linear stability of two inviscid fluids, in the presence of gravity, sheared 
past each other and separated by an flexible plate. Conditions for exponential growth 
of velocity perturbations are found as functions of the flexural rigidity of the plate and 
the shear rate. This Kelvin-Helmholtz instability is then analysed in the presence of 
plates with spatially periodic (with period a) flexural rigidity arising from, for example, 
a periodic material variation. The eigenvalues of this periodic system are computed us- 
ing Bloch's Theorem (Floquet Theory) that imposes specific Fourier decompositions of 
the velocity potential and plate deformations. We derive the nonhermitian matrix whose 
eigenvalues determine the dispersion relation. Our dispersion relation shows that plate 
periodicity generally destabilises the flow, compared to a uniform plate with the same 
mean flexural rigidity. However, enhanced destabilisation and stabilization can occur for 
disturbances with wavelengths near an even multiple of the plate periodicity. The sensi- 
tivity of flows with such wavelengths arises from the nonpropagating, "Bragg reflected" 
modes coupled to the plate periodicity through the boundary condition at the plate. 



1. Introduction 

When two immiscible fluids are sheared relative to each other, the flow, and the nor- 
mally flat interface between them, may become unstable. This classical Kelvin-Helmholtz 
instability ( |Chandrasekar (1961)| [Drazin fc Reid (2004)| [Whitham (1999)| ) arises from 
perturbations of the fluid velocity that grows into a vortex sheet due to local Bernoulli 
variations in pressure. The presence of gravity also tunes this instability, depending on the 
mass density difference and shear rate between the two fluids. If the denser fluid is in lower 
gravitational potential, there is a critical shear rate beyond which small perturbations 
grow exponentially. This instability first arises in modes of the fluid velocity perturbation 
that have wavevectors greater than a critical wavevector. All velocity perturbations with 
sufficiently small wavelength become unstable as long as the shear rate exceeds a criti- 
cal value. However, very small wavelength perturbations can be stabilised if additional 
restoring forces from, say, an interfacial surface tension (Whitham (1999) I, or a separat- 
ing elastic plate are included. The Kelvin-Helmholtz instability, and its generalisations, 
have been thoroughly studied, for example, in the context of knotting in astrophysical jets 



(Blanford & Rees (1974) I, cloud top entrainment instabilities (Reiss & Corona (1977) I, 
magneto-hydrodynamic instabilities ( Chandrasekar (1961) and Frank et al. (1996) I, and 
the flutter of panels | ,Miles (1956)| ). 

In the context of potential applications involving fiuid-structure interactions, we con- 
sider the Kelvin-Helmholtz instability in the presence of an infinitesimally thin bendable 
plate separating the two inviscid flows, as shown in Figure We assume the plate is 
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Figure 1. Discontinuous shear flow separated by a periodically flexible plate. 



thin and inextensible, such that only bending modes can be excited. Since elastic restor- 
ing forces stabilize short wavelength perturbations, instabilities arise only for perturba- 
tions that have wavevectors within a certain window of values. Furthermore, we con- 
sider an elastic membrane that has spatially periodic properties, such as mass density or 
bending rigidity, which might arise, for example, in an engineered composite structure. 
By itself (in the absence of fluid), a periodic plate supports stable bending waves that 
propagate except at band gaps, i.e. frequencies corresponding to waves that have half- 



integer wavelengths in one period of the structural variation ( Ashcroft & Mermin (1976) 
and Joannopolous, Meade, & Winn (1995)1. Waves near these wavelengths have van- 



ishing group velocities and do not propagate in an infinite periodic medium due to 
their multiple "Bragg" reflections and deconstructive interference. Properties of wave 
propagation through periodic media have been computed in numerous physical set- 
tings, including acoustic waves moving through materials with periodic sound speeds 



(Sigalas & Economou (1992) I, electromagnetic propagation through periodic dielectrics 



( Joannopolous, Meade, fc Winn (1995)|), and water wave pr opagation over periodic struc- 
tures dChen et al. (2004)||Chou (1998)||Hu fc Chan (2005)||Hire^ al. (2003)||Naciri fc Mei (1988)| 



Peter, Meylan, fc Linton (2006)[ [Porter fc Evans (1996)| ) 



In this paper, we combine the periodic stabilizing membrane with an otherwise unsta- 
ble shear flow and study the complex dispersion relation. A velocity perturbation with 
any wavevector, convected by the base laminar flow forms a travelling wave. This trav- 
elling wave, when coupled to the intrinsic periodicity of the elastic plate may be Bragg 
reflected and not able to propagate within the region of periodicity, altering its stability 
characteristics. In order to explore this heuristic argument further, we formulate in the 
next section the linear problem in terms of a vorticity-free velocity perturbation and a 
small deformation of the separating elastic membrane. We review, in turn, a uniform 
plate subject to destabilizing flow, and a free periodic plate that supports stable waves. 
The fully coupled problem, a periodic plate in destabilising flow, is then presented in 
the Results section. We show that the stability properties are affected most when the 
wavelengths of the most unstable modes are close to the periodicity of the plate flexibil- 
ity. Implications and extensions to related flow systems are discussed the Discussion and 
Conclusions. 
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2. Linear Model 

Consider the fluid-structure interaction depicted in Figure^ Two incompressible fluids 
flowing uniformly in two different directions, U+ and U_, in the x — y plane are separated 
by a negligibly thin elastic plate. Gravitational acceleration is in the — z direction. Since 
viscosity will be neglected, we will consider only potential flows. The total velocity in 
each region (±) is defined by the uniform flow U± plus a small irrotational perturbation 
v± = ^ip±. Upon imposing incompressibility, the upper and lower velocity potentials, 
ip±, satisfy Laplace's equation: 



(Vi+9>±(r,z,t) = 0, Vl = dl + dl (2.1) 

Since dzip\z=a ~ is the vertical velocity of the surface, the kinematic conditions at 
the plate-separated interface are 



dth{r,t)+V± ■ \/±h{r,t) - lim 9^(^±(r, ±e, i) = 0, (2.2) 

where h{r, t) is the normal displacement of the plate. The linearized normal stress balance 
at the interface, including the normal forces arising from the flexible plate, is 



[pip + plJ-W±ip + pgh{r, <)]+ = m(r)/i(r, t) + Wl{D{r)Wl)h{r, t), (2.3) 

where [X]^ = Ihn^^o X{z — s) — X{z — — e). In Eq. 12.31 m(r) and D{r) are the 
mass per unit area and the flexural rigidity of the plate, respectively. Since h(r, t) and 
</5(r, z, t) are defined only up to a constant, we need only consider the nonzero wavevector 
decompositions of the plate displacement and velocity potential, 



Mr,i) = E / ^(-)e---*^ and ^±(r,.,i) = J] / ^±(a.)e— ^l^-'^'^^. 

(2.4) 

In Eas. 12.41 we have also assumed a continuous time-harmonic decomposition. The rep- 
resentation of (^±(r, 2;,t) automatically satisfies Eg. 12.11 and the boundary conditions at 
z ±oo. 

The periodic plate mass density and rigidity functions that obey m(r) = m(r + R) 
and D{y) = D{r + R), where R = niaiei -I- 7120262 (ni, 712 integers) is any basis vector 
describing the periodicity. Here, ai, 02 are the periodicities of ?7i(r) and D{r) in the ei, 62 
directions. The periodic functions m(r) and D(r) can thus be decomposed into a Fourier 
sum over reciprocal lattice vectors G, 

m(r) = ^ m(G)e'*^ '" and i?(r) = ^ i?(G)e*° '". (2.5) 

G G 

The reciprocal lattice vectors G = tJiGi -1- n2G2 are defined by the reciprocal lattice 
basis vectors 



27r ea X z 27r z x ei 

Gi = and G2 = -. (2.6) 

oi ei • (e2 X z) a2 ei • (62 x zj 

Finally, we assume a large plate with dimensions Li x L2 that fit an integer multiple 
of the period (Li = Niai and L2 — N2a2) and impose periodic boundary conditions 
on the entire system: m{Y -f NiGiBi) = m{r) and D(r -|- NiaiBi) = D{r) for i = 1,2. 
Under this construction, the wavevector sum in the decompositions 12.41 are taken over 
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q = X)n^o("'«/-^')*^«- the Li, L2 — > cxd limit, sums over q become the Cauchy principle 
value of the corresponding Riemann integral. 

After substitution of Eas. ini into the kinematic conditions l2.2l we find the relationship 



(2.7) 



Upon substitution of Eas. lT^ 12.51 and 12. 71 into Eg. 12.31 and exploiting the orthogonality 
of the e^'^ '' basis functions, 

[ip+ + p-)u;^ - 2w(p+U+ + p_U_)-q + p+(q-U+)2 + p_(q.U_)2 + {p+ - p_)g|q|] 



- -^1q| Eg' "^(g')^-g' + Eg' Iql1q - GfDiG')h^^G'- 

(2.8) 

Equation 12.81 shows that each coefficient /iq is linked with Ni x N2 other coefficients 
/iq_G'- This problem can be solved numerically by truncating /iq_G ~ for large |G| 
where to(G) and D{G) are small. Alternatively, Eg. 12.81 can be organized according to 
the "reduced zone scheme," by defining k = q — G, where G is such that k is restricted 
in the first Brillouin zone defined by |k • e^j < |Gi • 6^1/2 (Ashcroft & Mermin (1976) I. 
In this representation, we need only consider the A'^i x A^2 eigenvalues, each with index 
k inside the first Brillouin zone. The problem is thus cast into a quadratic eigenvalue 
problem defined by 



E 

G' 



[Ak(G, G>2 + B^{G, G')lu + Ck(G, G')] /ik-G'O 



UJ 



0, 



where 



Au{G,G' 
Sk(G,G') = -2 
Ck(G,G') 



^.o,G^+^(G'-G) 



U+-(k-G) U_-(k-G) 



|k-G| 



(U+-(k-G))^ 
^+ Ik-GI 



" |k-G| ^ 
(U_.(k-G))2 



^^G^G' 



+ (P+ - P-)9 



|k-G| 
|k- G|2|k- Gf D(G' - G 



3G,G' 



(2.9) 



(2.10) 



Eouation l2.9l is a matrix equation for hu^G{uj) that has solutions for only certain u;(k). 
The full dispersion relation can be recovered by assigning the appropriate translation G 
in the eigenvalues a;(k — G) to tt'(q) for wavevectors q outside the first Brillouin zone. 
The system becomes unstable when w(q) acquires an imaginary component. 

In anticipation of a length scale a associated with the periodicity in the plate properties, 
we henceforth nondimensionalise all parameters and eliminate p- , g and a according to 



r r/a, G aG, k ak, uj 
D 



LU, U± 



ag 



D 



m 



(2.11) 



a'^p-g' 



m 



, and hk(G) ^ hk(G)/a. 
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Figure 2. (a) The curve relating the critical value of the shear (7* to the constant flexural 
rigidity Do, for various values of the uniform plate mass density mo. For velocities greater than 
[/, , the flow becomes unstable. For different plate mass densities, the sets of D, , )7, that give rise 
to an initial instability at wavevectors — tt and = 27r are shown by the black arcs spanning 
the curves, (b) The dispersion relation for a uniform plate just after the onset of instability. 



2.1. A Uniform Plate in the Presence of Flow 
First consider the standard Kelvin-Helmholtz instability in the presence of a uniform 



elastic plate, where m = mo and -D(r) = Dq are constants. Miles (2001) treated a 
similar problem of a thin boundary layer over an elastic plate. Since D{G) = Do^g.o and 
m(G) = rnQScfl, Eg. 12.91 is diagonal and is satisfied when 



r(gV±(q) = (7U+ + U_).q±yA(^, (2.12) 
where r(g) = 7 + 1 + moq, and the discriminant 

A(q) = r{q){Doq' + (1 - j)q) - - lJ^)-qf - moq{^{lJ+-qf + {U^-qf). (2.13) 

The eigenvalues lu± become complex when A(q) becomes negative, leading to exponen- 
tially growing perturbations. In the limit where the influence of the plate is negligible 
{mo, Do « 0) the dispersion relation reduces to that of the standard Kelvin-Helmholtz 
instability. The inclusion of membrane stiffness (Do > 0) tends to stabilise small wave- 
length modes. For any given mo ^ 0, the pairs of critical Do and U (denoted D^ and ?7*) 
that give the onset of instability at a chosen wavevector can be expressed explicitly in 
two flow configurations where U± and the q corresponding to the most unstable velocity 
perturbations are all collinear. When U± are collinear, the onset of instability arises 
when A(q) = is solved by two critical repeated, real roots, g*, in addition to the root 
at (? = and two complex conjugate roots. For a static configuration of the lower fiuid 
(U_ = in the reference frame of the plate), we find 



(1 + mpq^y^ - 7^(1 + 2mog^) - 'ym'^ql 

g4(3(l + 7 + "^o'7*) + 2(7 + 3)mog*) ^ ^ 

(2.14) 

4{1 - + J + moq,f 



D, 



7g*(3(l + J + mggj) + 2(7 + 3)mog,) ' 
Alternatively, for a chosen value of Do, the critical t/* and the corresponding initially 
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unstable wavevector can be found implicitly. The critical velocities as a function 
of Do and mo are shown in Fig. El^a). Note that the critical velocity J7* decreases as 
the plate mass density mo increases. The limit mo — > cx) is different from the p± — > 
limit treated in the next section. Here, the plate flexibility remains constant as its mass 
density increases. Although the magnitude of uj±{q) decreases, the discriminant (Eq. 
I2.13|l becomes negative at smaller [/» due to the increased inertial forces of he plate back 
on the fluid. 

The two branches of the full dispersion relation ijo±{q) are plotted in Fig. Olb). We 
have chosen parameters (_D»(7r) and defined by Eas. 12.1^ such that the onset of 
instability arises at « tt. For U = l/v^ < U^{q^ — n) = 0.756320, uj±{q) are real, 
as shown by the dashed curves. We also plot (solid curves) the dispersion relation for 
U ~ 0.76f577 > = tt) = 0.75632, whereupon a window of instability opens near 

= TT. The magnitude of the complex parts of uj±{q) are indicated by the height of the 
shaded region. 

In the second flow configuration, we assume a vanishing net momentum in the reference 
frame of the plate: U_ = — 7U-)_. In this case, we find independent of mo, 

1-7 9 4(1-7) , . 

D*^^ and U! = ^ f ^ ' , 2.15 
3qt 37(7+1)9, 

and the symmetry uj+{q) = —uj-{q). For Dq < and/or J7 > ?7,, as defined by Eq. 
12. 141 or Eg. 12. 151 ui{qfv q^,) becomes complex and instability arises for a small window of 
wavevectors near . 

2.2. Dispersion Relation for a Free Periodic Plate 

Now consider the limit of an isolated (but periodically structured) plate where the inertia 
of the bounding fluid is negligible. If the plate were uniform, the dispersion relation in the 
limit p± ^ is simply mow^ = D^q'^ . For a periodically structured plate (m(G 0) ^ 
0, D(G 7^ 0) 7^ 0), the dispersion relation can be found from Eg. 12. 91 with the simplified 
matrices A{G - G') = m(G - G'), S(G - G') = 0, and C(G - G') = -|k + Gpjk + 
G'pZ)(G — G'). Henceforth, we assume the periodicities in the plate to be sinusoidal in 
the x-direction and take the forms 



m(r) = mo + 2mi cos 27rx; D(r) = _Do + 2Z?i cos 27rx, (2-16) 

Although the mean mass density and bending rigidity are unchanged from mo and Do, 
respectively, the variations mi and/or Di will change the dispersion relation. Specifically, 
periodic variations will break the dispersion relation at wavevectors corresponding to 
waves that cannot propagate through an infinite, periodic medium. At each of these 
wavevectors (the "Bragg planes"), the waves (half-integers of which fit into a period) 
become standing. Corresponding to these nonpropagating waves are gaps in the frequency 
(band gaps, or stop bands) within which an infinitely periodic material cannot be excited. 

Figures ^^,h show the dispersion relation of an isolated plate with periodicities in the 
bending rigidity (a), and mass density (b). The frequency gaps increase as the contrast 
rui or Di is increased. From consideration of the eigenvectors, the standing modes near 
the band gaps can be shown to have either nodes or antinodes over the high flexural 
rigidity (or high mass density) regions, depending on which side of the band gap the 
wavevector q lies. These behaviour generally arise in wave propagation through periodic 
media. When the plate is coupled to flow, it will be the near equal-wavelength modes 
immediately straddling the band gaps, with different frequencies and ~ 90° out-of-phase, 
that will qualitatively affect the stability near the band gaps. 




q q 

Figure 3. (a) The normalized dispersion relation ^J'moJTJo Lu±{q) for mi = 0, D\/Do = 0.1 
(dashed curve) and Di/Dq = 0.25 (solid curve), (b) A blow-up of the dispersion relation for 
Di = 0, mi/mo = 0.1 (dashed curve) and mi/mo = 0.25 (solid curve). 



3. Results for the Coupled Problem 

We now combine the known results of the previous sections and consider the fully 
coupled problem where fluid instabilities and Bragg reflection interact. For simplicity, 
we consider the plate periodicity wavevectors (G = 27rn/ax) and the uniform flows 
U± ~ U±x to be aligned. In the full problem, the dispersion relation uj{q) may be 
both complex (signalling regions of instability) and discontinuous at band gaps where 
{diu/dq) = 0. The full dispersion relation is found from solving the quadratic eigenvalue 
problem given by Eq. 12.91 This problem can be expressed in linear eigenvalue form, 



(M - t^I) • J7 = 0, (3.1) 

with T] = {hu,ojhu)^ , and 



M = . (3.2) 

The dispersion relation a;(k) is found from the eigenvalues of the nonhermitian matrix 
M. For the parameters explored, we find numerical convergence of the lowest handful of 
eigenvalues with 40 or fewer modes G, G'. Therefore, we truncate the system at approx- 
imately 60 modes (where M is a 120 x 120 matrix). This is more than sufficient to obtain 
numerical accuracy for the lower eigenvalues at all relevant wavevectors. The eigenvalues 
of M are sorted from the numerical solutions and replotted in the extended zone scheme 
where the wavevector ^ |q| < cx). When unfolding the eigenvalues w(k), we make use 
of the fact the functions a;(0 ^ |q| < oo) in the normal extended zone scheme come in 
pairs w±(q) and have the symmetry property ti'±(— |q|) = — Wzpdqj). This property, along 
with the consideration of the corresponding eigenvectors at each wavevector k allows us 
to reconstruct U!±{q) from a;±(fc). 

In Figures 0fa-c) we show the effects of increasing the strength of the periodic rigidity 
Di/Dq for the case U_ = 0. In Fig.QJa) we show a system with an incipient instability at 
g « g* = 7r/2. As Di is increased (Figs.0{b-c)) the bubble of instability grows modestly, 
while gaps in the dispersion relation also appear. The destabilising effect of rigidities 




q q q 

Figure 4. (a-c) The complex dispersion relation for a system just above the instabil- 
ity threshold. The instability is first triggered at q, — n/2 (which arises for parameters 
7 = 1/2, mo = 0.2, Do ^ D, = 0.0298349, and > = 1.19809). As Di increases, the 
instability bubble (Im{aj±}) increases in size as more wavevectors become unstable, (d-f) Dis- 
persion relation for instabilities that first arise at q, = 3.14 < tt. As Di is increased, the modes 
with q < n are destabilised, while those with q > tv are stabilised, (g-i) Increasing Di when the 
incipient instability starts at q, — 3.25 > tt destabilises modes with q > tv and stabilises those 
with q < n. 



with periodicity a 3> li^ j q-i, can be understood in terms of a slight effective softening of 
the plate to bending modes with wavelengths greater than af. 

Note that the gaps here do not appear at the Brillouin zone edge q = tt defined 
by a system's static periodic structure as in the free plate problem (Figs. and in 
most other settings. Rather, they appear symmetrically around the band edge due to a 
doppler shift arising from the nonvanishing net momentum of the system. Gaps occurring 
away from the zone boundaries have been shown by Bermel & Warner (2002) to arise in 
electromagnetic propagation through deformed cholesteric elastomers. 

Now consider an incipient instability that arises at higher wavevectors, near the first 
band gap at q ~ ir. Assume for simplicity, the zero net momentum (U_ = — 7U+) 
case where Bk(G,G') = in Eg. 12.91 and uj±{q) ~ —u!^{q). Figures 01 d-f ) show the 
effects of increasing Di when the initial instability starts at <?, — 3.14 < tt. As the 
plate heterogeneity is increased, modes with wavevector q < tt are destabilised, while 
those with q > tt are stabilised. Conversely, if the incipient instability arises at g, > tt, 

f The effective long wavelength bending rigidity of a plate made from al- 
ternating equal strips of two materials with individual rigidities Do ± AD is 
D,ff = 2/(l/(Do + AD) + l/(Do - AD)) ^ Do - {ADf/Do 
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as in Figs. 2fg-i) (where = 3.25), increasing Di will destabilise modes with q > tt, 
while stabilising those with g < tt. The discontinuous change in the eigenvalues across 
the Bragg plane q = n occurs first in the imaginary component and the growth rate 
Im{uj{q)} is discontinuous. As the contrast Di/Dq is further increased, Ke{u!{q)} also 
becomes discontinuous. This behavior can be explicitly seen with perturbation analysis 
of Eg. 13.21 bv truncating the matrix and finding the four eigenvalues as g — + tt. 



4. Discussion and Conclusions 

We have analysed the properties of the instabilities of two fluids uniformly flowing 
past each other, separated by an elastic plate with periodic bending rigidity and/or 
mass density. Our generalization of the Kelvin-Helmholtz instability can be quantified 
by finding the eigenvalues and eigenvectors of a nonhermitian matrix (Eq. 12. 9|) . The 
imaginary components of the eigenvalues uj{q) determine the rate of unstable growth 
for eigenmodes of wavevector q. We find that an imposed plate periodicity generally 
destabilises the flow but that the destabilisation mainly occurs in the band of wavevectors 
on the side of the initially most unstable mode. Modes with wavevectors on the other side 
of the band gap from the most unstable wavevector g* are stabilised as Di is increased. 

The destabilisation/stabilisation that occurs on either side of the band gap can be 
understood in terms of ~ 90° out-of-phase standing modes with wavevectors q = {2n — 
l)7r±£ immediately straddling the n*'* band gap. For example, if > (27i— l)7r, increasing 
Di destabilises the waves with q Ki q^, since their antinodes occur predominantly over 
regions of lower rigidity. Modes with q < {2n — l)7r, on the other side of the n*'' band 
gap, can be significantly stabilised since they are approximately 90° out-of-phase, and 
their antinodes sample a stiffer plate. The same argument applies for < (2n — l)7r. 

Our results suggest possible stability control strategies. By adding restoring forces 
(such as bending rigidity) with periodicity near the wavelength of the most unstable 
mode, stability can be enhanced or diminished at wavelengths either less than or greater 
than the band gap wavelength. Stabilisation of specific modes against Kelvin-Helmholtz 
instabilities can be achieved, but at the expense of destabilisation of other nearby modes. 

Extensions of our approach to stability analysis in other fluid-structure systems should 



be possible. For example, the Orr-Sommerfeld problem lOrszag (1971) I in nonuniform 
Poiseuille flow through periodically elastic pipes can be similarly solved. A physical re- 
alization of flow through elastically periodic pipes may have arisen in the experiments of 



Kramer (1960) in which damping elements were periodically embedded along an elastic 



pipe. The Kelvin-Helmholtz type instability has been only briefly mentioned in this con- 



text by Benjamin (1960) and Landahl (1961) and only in uniform pipes. Moreover, the 



effects of periodicity on other classes of instabilities {e.g. Tollmien-Schlichting) have not 
been explored. It would also be interesting to consider the effects of dissipation, in both 
the fluid flow and the viscoelastic plate, on the stability properties of a periodic system. 

The author was supported by the US National Science Foundation through grant 
DMS-0349195, and the US National Institutes of Heahh via grant K25AI058672. 
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